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Abstract 

We calculate the Regge trajectories of the subleading BFKL singularities and eigenfunc- 
tions for the running BFKL pomeron in the color dipole representation. We obtain a viable 
BFKL-Regge expansion of the proton structure function F 2p (x, Q 2 ) in terms of several right- 
most BFKL singularities. We find large subleading contributions to F 2p (x, Q 2 ) in the HERA 
kinematical region which explains the lack of a predictive power of GLDAP-extrapolations 
of F 2p (x, Q 2 ) to a domain of small x. We point out the relation of our early finding of the 
precocious BFKL asymptotics to the nodal structure of subleading BFKL eigenfunctions. 
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Extrapolation of the proton structure function from the accessible region of x and Q 2 
to small x remains the topical issue ever since the first SLAC experiments on deep inelastic 
scattering (DIS). The GLDAP evolution has been a standard instrument in these studies. 
It was soon realized that GLDAP extrapolations are not unique, and equally good fits to 
the large- x-data did invariably, and wildly, diverge with each other and new experimental 
data, when extrapolated beyond the studied range of x. For instructive comparison of 
extrapolations of the pre-HERA fits to the proton structure function with the new data from 
HERA down to x ~ 10~ 5 see 0. The lack of predictive power at small x is not surprising: 
the GLDAP evolution describes the future for large Q 2 starting with assumptions on the 
parton densities at past, Qq. One needs to assume the whole function of x and different 
past defines a different future. 

At very small x, the assumptions of the GLDAP evolution break down and it is su- 
perseded by the log(l/x) BFKL evolution ||. The BFKL evolution is meant to predict 
the future at small x from the past defined as a function of Q 2 at a sufficiently small xq. 
At an asymptotically small x the BFKL prediction, F 2 (x,Q 2 ) = Fjp(Q 2 )x~ Alp , is unique 
modulo to the total normalization factor. In the BFKL approach, the x-dependence of the 
structure function at moderately small x depends on the spectrum of eigenvalues and on the 
eigenfunctions of the BFKL operator. The purpose of the present studyf] is an evaluation 
of contributions from subleading BFKL singularities for the running BFKL pomeron in the 
color dipole formulation || |, 

The BFKL equation for the interaction cross section section cr(£, r) of the colour dipole 
r with the target reads (hereafter £ = log(l/x)) <9cr(£, r)/<9£ — K ® <r(£,r). Here the 
kernel K, is related to the wave function squared of the color-singlet qqg state with the 
Weizsacker- Williams soft gluon. The gluon fields are calculated with the running QCD 
coupling and perturbative gluons are infrared regularized imposing a finite propagation 
radius R c ~ 0.2 — 0.3 fm. The BFKL pole IP n corresponds to the Regge-behaving solution 
of the BFKL equation, 

(Xq \ An 
—) , (1) 

where the eigenfunction a n (r) and the eigenvalue (intercept) A n are determined from 

K ® a n = A n a n (r) . (2) 
The behavior of eigenfunctions at r — > and/or ats(r) — *> has been found in || [7|: 

4 

, 7 n A n = - . (3) 

Useful clues come from the experience with solutions of the eigenvalue problem for the 
Schrodinger equation. First, the leading eigenfunction cr (r), found numerically in |5|, ||, [7], 
is node free. Consequently, the subleading solutions with eigenfunctions A n < A must 
have nodes. Second, the asymptotics ([3]) must hold for all eigenfunctions in the region of r 
beyond all nodes. Third, because of the infrared regularization cr n (r) saturate at r > R c . 
This suggests that for eigenfunctions with n nodes one must try a n {r) = V n (z)<j (r) where 

1 The preliminary results have been reported at the DIS'97 Workshop 
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V n {z) is a polynomial of a suitably chosen variable z ~ [l/as(Q 2 )]' y ■ Then we apply the 
variational principle and minimize the functional 



dr 



JC <g) aJr) - A n a r . 



(4) 



in the space of trial functions with the polynomial prefactor V n {z). The success of such an 
unorthodox variational principle for excited states depends on the guessed trial function, 
examples of astonishingly successful applications can be found in ||. 

With the so obtained eigenfunctions a n (r) we can calculate the contributions F^ L \Q 2 ) 
to the proton structure function using the color dipole factorization || and can perform the 
BFKL-Regge expansion of the proton structure function 

x Q \ An 



F 2p (x,Q 2 ) = Y,A n F^(Q 2 ) p . (5) 

Our principal findings on solutions of the BFKL equation are as follows. 
The running BFKL equation gives rise to a series of moving poles in the complex j-plane. 
The intercepts A n (Figure 1) very closely, to better than 10%, follow the law 



suggested earlier by Lipatov from the quasiclassical consideration [TO]] .The found eigen- 



functions plotted as S n (r) = a n (r)/r (Figure 2) to a crude approximation are similar to 
Lipatov's quasi-classical eigenfunctions, which are £ n (r) ~ cos[0(r)] for n ^ 1 fill. For 



a related numerical analysis of the running BFKL equation see [JTTJ] . Within our specific 
infrared regularization the node of 0"i(r) is located at r ~ 0.05 — 0.06 fm. With growing n, 
the location of the first node moves to a somewhat larger r, and the first nodes accumulate 
at r ~ 0.1 - 0.15fm. 

The slope a' n of the Regge trajectory for the pole IP n can be found using the technique 
of Ref.p], HH. We find a' n ~ 0.07 - 0.06 GeV~ 2 with weak dependence on n. 

The structure functions F^ 1 \q 2 ) for the IP„ poles are shown in Figure 3. At large 
Q 2 , far beyond the nodal region, F£(x, Q 2 ) oc (x /x) An [l/as(Q 2 )] 4 ^ 3A '\ Since the relevant 
variable is a power of the inverse gauge coupling, which varies with Q 2 very slowly, the 
nodes of F^iQ 2 ) are spaced by 2-3 orders of magnitude in the Q 2 -scale and only the first 
two of them are in the accessible range of Q 2 . The first nodes of F^iQ 2 ) are located at 
Q 2 ~ 20 — 60GeV 2 . Hence, only the leading structure function contributes significantly in 
this region. This explains the precocious BFKL asymptotics for Q 2 ~ 60 GeV 2 found earlier 



from the numerical solution of the color dipole running BFKL equation (Tj . 

An interesting finding is that the Born approximation for the dipole cross section, <Jb(t), 
gives a very good quantitative description of the proton structure function at x ~ 0.03 jTJJ 
The r-dependence of cr B (r) is quite different from that of the leading eigenfunction o"o(r), 
and expansion of cr s (r) in BFKL eigenfunctions shows that the contribution of subleading 
terms with n > makes up to ~ 60% of <Jb{t)- Consequently, the Q 2 -dependence of the 
proton structure function at x = xq, and the subsequent x-dependence towards smaller x, 
is controlled to a large extent by the subleading contributions F^iQ 2 ). 

At small x only Q 2 < 10 3 GeV 2 are accessible. In this range the structure functions with 
n > 3 are hardly distinguishable. Besides, the splitting of intercepts with n > 3 is much 
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smaller than for n — 0, 1, 2. Hence, the Regge expansion (|5]) can be truncated at n = 2 and 
F^\Q 2 ) comprises contributions of all poles with n > 2. 

The BFKL equation allows one to determine the intercepts and structure functions 
F^ l \Q 2 ). In the expansion we put A Q = 1. Then the overall normalization of eigenfunctions, 
for instance, a n (r ^> R c ) = 0.89 mb, and A\ = 0.39 and A 2 = 0.33 are two adjustable 
parameters which are fixed from the boundary condition at x = xq- With the proper 
account of the valence [15] and non-perturbative Jl4| corrections to (|) we arrive at the 



three-pole-approximation which provides a viable description of the experimental data [[16 
in a wide kinematical range (Figure 4). Notice that in the pre- nodal region of Q 2 < 20 GeV 2 
the leading and subleading structure functions are very close in shape and the experimental 
data in such a limited range of Q 2 are absolutely insufficient for the determination of A n , 
which explains the failure of the early GLDAP fits. 
The effective pomeron intercept 

d\ogF 2p (x,Q 2 ) 

A e// — { ' ) 

O log X 

gives an idea of the role of the subleading singularities. The intercept A e // calculated with 
the experimental kinematic constraints is much smaller than A = 0.4 which is expected to 
dominate asymptotically. The agreement of our numerical estimates with the HI determi- 
nation [[nj (Figure 5) is quite satisfactory. 
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Figure captions 



Fig.l BFKL eigenvalues. 

Fig. 2 BFKL eigenfunctions. 

Fig. 3 Modulus of BFKL structure functions. 

Fig.4 Three-pole approximation vs. the data of HI, ZEUS and E665. 



Fig. 5 Effective intercept vs. HI determination |L7 
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